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Abstract 



> 

C^ . The early stages in the evolution of the gluons produced in the central 

region of a head-on high-energy heavy ion collision are studied. An equation 



(N 

m 



^ [ is given for the rate of change of transverse momentum into longitudinal 



O 



momentum where the longitudinal direction is along the collision axis. We 
0> ' are able to follow the system up to the time where equilibration seems to 

^ ■ be setting in, but we are unable to actually follow the system as it reaches 

^. equilibrium. 

Oh! 

• . : 1 Introduction 

>. 

^ ■ At very small values of x the gluon density in a proton's light-cone wave- 

j^ . function reaches saturation [|l|, ||, H, ^, that is for gluons having transverse 

momentum below the saturation momentum, Qs, there is a density on the 
order of ^ ^~ ' gluons per unit of transverse phase space (transverse coor- 
dinate times transverse momentum space). Parton saturation in the proton 
may already have been observed at HERA[^, |^. This regime of high density 
gluons is reached at more moderate values of x for a large nucleus. In a 
head-on collision of high energy heavy ions saturated gluons, those having 
transverse momentum at or below Qs, are freed at very early times after the 
valence quarks of the two ions have passed each other, say in the center of 
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mass frame. It is widely assumed that this gluon system then quickly reaches 
kinetic equilibrium and somewhat more slowly reaches equilibrium between 
gluons and quarks. 

There are Monte Carlo calculations which have followed the early stages 
of the evolution of a heavy ion collision and which give good evidence for 
equilibration]^, §]. There is also an interesting Monte Carlo calculation0 
which followed the time evolution of the gluon fields which start from the 
saturated distribution occurring in the McLerran-Venugopalan model. Our 
purpose here is somewhat different. The object is to give a rough analytic 
calculation of how gluons evolve in momentum and in space after being freed 
in a high energy heavy ion collision. We are able to follow a typical gluon 
through its many small angle scatterings up to the point where the interac- 
tions which should actually give equilibration begin to occur. At the time 
where equilibration actually begins to set in our technique breaks down, al- 
though we are able to say at what time and at what momentum, relative to 
the saturation momentum, equilibration starts to set in. 

The starting point of our discussion is the McLerran-Venugopalan |I0| 
model for the light-cone wavefunction of a heavy ion|jll|. This model is 
characterized by one parameter the saturation momentum, Qs, below which 
gluon densities reach their maximum value. We suppose that a finite fraction 
(our parameter c in (16)) of all gluons having transverse momentum less than 
or equal to Qs are freed in a head-on heavy ion collision [|1^, and that those 



gluons having small rapidity are freed in a time about equal to 1/Qs- Gluons 
having small rapidity can then scatter only with other gluons having small 
rapidity, because the rate of scattering is slow and higher rapidity gluons 
quickly separate from lower rapidity gluons. It is then straightforward to 
write an equation for the rate of change of the transverse momentum as zero 
rapidity gluons scatter. The equation is given in (49). In arriving at (49) 
we have also had to determine where the infrared cutoff lies for small angle 
scattering. We have done this by requiring that the soft gluon field, coming 
from the hard gluons we are dealing with, also be limited in magnitude to 
be no larger than 1/g. This determines a minimum m 
omentum transfer, £„, which is not quite the same as has been suggested 



earlier] 13]. The difference between i^ given in (44) and the screening length 
of Ref.l3, given in (45), changes the constants in (49), but not the form of 
the equation. 

What has been said for zero rapidity gluons also applies to gluons having 



rapidity y, by simply considering their interactions in a frame boosted by -y. 
Thus gluons of rapidity y only interaction with other gluons having rapidity 
y and (49) remains valid for the rate of change of transverse momentum if 
i^i + lnchy. 

At early times the rate of change of transverse momentum is much less 
than that corresponding to the one-dimensional expansion of a gas of gluons 
in equilibrium. However, because of the time dependence of the infrared 
cutoff the rate of loss of the square of the transverse momentum of gluons 



increases so that at a time, for zero rapidity gluons, on the order of ^ eV ''"^^ , 

and at a transverse momentum given by {^l/QsY ~ V ^f^' ^^^ expansion of 
our gluon gas more closely resembles that of a system in equilibrium. At this 
same time scatterings which change longitudinal momentum into transverse 
momentum, not included in (49), begin to become important. This is also 
the time at which gluon interactions no longer take place simply among 
gluons of the same rapidity so that our whole formalism is breaking down. 
Finally, this is the time at which the z— extent of the medium is becoming 
greater than the mean free path for gluon interactions which can change a 
gluon's momentum by an amount comparable to its momentum. It seems 
very likely that this is the time at which equilibration is setting in, although 
we are unable to follow the system further because of the limitations of our 
present formalism. It is important to note that the time when equilibration 
seems to be setting in is parametrically short compared with the radius of 
the initial colliding ions. We consider this evidence that equilibrium indeed 
occurs in very high energy relativistic heavy ion collisions while the system 
is still undergoing a one-dimensional expansion. 



2 Small-x wavefunctions and parton satura- 
tion 

In this section we shall review the parton description of the light-cone wave- 
function of a heavy ion. There is a very simple model, the McLerran- 



Venugopalan [jlOi model, which very nicely illustrates parton (gluon) satu- 
ration. Indeed, it has recently been claimed 0] that the main results of the 
McLerran-Venugopalan model should be general results in QCD. We shall 
briefly review the physical basis of this model and then comment on the 



general expectation of the reliability of the model. 

We begin, following McLerran and Venugopalan, by looking at the distri- 
bution of valence quarks in a high energy ion. The valence quarks are found 
in a Lorentz-contracted longitudinal disc of size A z = 2R ■ — where R is the 

O p 

nucleon size and p the momentum per nucleon of the ion. For purposes of 
our qualitative discussion we consider A z = so that one can imagine the 
valence quarks having a two-dimensional number density of 



ng{b) = 6pVW^^ (1) 

valence quarks per unit area where p is the normal nuclear density while b is 
the impact parameter measured from the center of the nucleus in a direction 
perpendicular to the direction of motion of the nucleus. 

The valence quark density given in (1) is the source of soft gluons corre- 
sponding to the Weizsacker- Williams field of rig^ [1^. The color charge at a 
given impact parameter comes from a random addition of the color charges 
of each of the valence quarks at that impact parameter. A single quark gives 
^^ inQ'^/n'^ gluons at scale Q^ and per unit rapidity so that one expects 
the number density of gluons per unit area and per unit rapidity to be 

dxGA{x,Q ) ru\^^P B na I 'i M\ 
^ = ^gW—^ ^^ Q //^ (2) 

leading to a gluon distribution in the nucleus of 

xGa{x, Q') = SA'^-^in Q^fi^ = AxG{x, Q^). (3) 

IT 

(We note that for a proton (3) is not too bad when 10^^ < x < 10^^ and 
Q^ = 5 — lOGeV^.) The unitegrated gluon distribution 

dxGAix,P) /i^-^ctCp ,,s 

has the interpretation, in light-cone gauge, of the number of gluons per unit 
rapidity per unit of transverse phase space (impact parameter space times 
transverse momentum space). Of course (2)- (4) are correct only in the low- 
density limit. In the high density regime the distribution of gluons in the 



wavefunction of the nucleus in the McLerran- Venugopalan model isg, |11 



dxG N —1 f d Xi ~x'^0'^/A\ il-x /r\ 



d%d'^i± ATT^aNc J x} 
where the saturation momentum, Qs, is given by 

Ql=^-^^4W^- pxG{x^^). (6) 

The large (^ hmit given in (4) comes from the logarithmic singularity of 
xG when evaluating (5). In the region (?i_/Ql << 1 one can neglect the x^ 
dependence of Ql and write 

dxG Ar2-l foo ^ ^_^^,^^^, 



with the scale of xG in (6) taken as t- = Q^. From (7) one finds, when 

il/Q's « 1, 

'"'' ""^-^ in Ql/il (8) 



which exhibits saturation, the logarithm naturally present in the momentum 
space expression. In Ref.4 a result essentially identical to (8) has been derived 
using BFKL dynamics. We expect the essential features of (7), and (8), to 
be general though it is difficult to reliably set the overall constant in (8). 

Despite our notation it should be emphasized that (8) refers to the dis- 
tribution of gluons in a light-cone wavefunction. It has no relation to the 
gluon distribution, or its derivatives, as determined by an operator product 
expansion analysis of deep inelastic lepton-nucleus scattering. The result 
in (8) does, however, give the transverse momentum distribution of gluons 
produced in a hard scattering off a large nucleus. 

Of course, the whole approach discussed here only makes sense if the 
parameter Q^ is in the perturbative regime. For RHIC energies (6) gives 
Qs ~ IGeV while for LHC energies Qs should be 2 — 3GeV. Thus, our 
perturbative approach is certainly marginal for central RHIC collisions while 
it should be more applicable for heavy ion collisions at the LHC. Of course 
the same phenomenon of saturation also occurs at very small x values for a 
proton and, indeed, gluon saturation may have been observed at HERA. 



3 Toward equilibration in the central region 
of a head-on heavy ion collision. 

In this section we shall discuss, in a semiquantitative manner, the evolution 
from very early times up to the time where equilibrating interactions begin 
to occur for the central region of a zero impact parmeter heavy ion collision. 
There is good evidence from Monte Carlo simulations that equilibration does 
indeed occur in collisions having kinematics similar to those which will take 
place at RHIC and at the LHC[^, If. Our purpose here is of a more theoretical 
nature, to try and understand in analytic detail the early stages of evolution 
of a dense gluon system and to see that equilibrating interactions do begin 
well before the system falls apart. 

3.1 Which gluons are freed and when? 

At the time of collision of the two ions the gluons in the wavefunctions of 
the ions are distributed according to (5). For (^ IQ\ << 1 the distribution is 
given by (8), a distribution which is, except for a logarithm, given by phase 
space. Since there are few gluons in the wavefunctions having (-^/Ql << 1 
we shall simply ignore them completely in determining the initial distribution 
of gluons. When (!-\/Ql >> 1 the distribution of gluons in the wavefunctions 
is given by (2) and (4). These gluons are additive in the various nucleons 
in the nucleus so that one may write the two gluon (jet) production cross 
section as 



da f 2, j2u dxiGAixif) dx2GA{x2,f) da 



dyidy2d'^i± d% J d%i d'^b2 dt 

(9) 
where xGa is given by (2) and where the gluon-gluon elastic scattering cross 

section is 

da aNc.2 ^vr^ ut us st 

~di ^ ^^T^ (iv2- 1)52^3 - ^ - ^ - ^) (10) 

with t = — £^. The main interest of (9) and (10) for us here is that for 
^i/Qs >> 1 the cross section in (9) scales like ij_'^. Thus one gets more 
particles, and more transverse energy, by lowering i±. Thus the dominant 
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region for gluon production, and for the production of transverse energy, in a 
very high energy heavy ion colhsion is the region where the produced gluons 
have transverse momentum on the order of (5s;[]T2| ^^^ where the gluons 



come from the interaction of gluons having transverse momenta on the order 
of Qs in the wavefunctions of the colliding ions. 

Now that it is clear that the main contribution to the initial transverse 
energy density in a high energy heavy ion collision comes from gluons having 
transverse momentum on the order of Qs, the saturation momentum, we 
need to ask at what time are these gluons produced. Suppose we focus on a 
unit of rapidity centered about Yq, —^ + Yq < y < ^ + Yq with quanta in this 
region being right-movers in the frame of choice. Then define P by 

P = Qs sink Yo. (11) 

P is the typical z-component of the momentum of the gluons we are con- 
sidering. The time over which the gluons in our unit of rapidity are freed 
is 

f is the time over which gluons in differing units of rapidity physically sepa- 
rate with those we are considering. It is also the time over which our gluons 
physically separate from the valence quarks. We may estimate the time of 
this latter separation by neglecting the z-extent of the right-moving valence 
quarks. The z-extent of the gluons in — ^ + Fq << 1/ < ^ + ^o is 

A^ = i (13) 



P 



while their longitudinal velocity is 



- 1 - -^. (14) 



p2 + Q2 2P2 

The time which it takes the right-moving valence quarks to separate from 
the gluons in our unit of rapidity is 

At = ^=^. (15) 



Using (13) and (14) we see that Ar = f. 

Thus the gluons on which we are focusing are freed during a time given 
by (12). We have described these gluons as being virtual gluons in the light- 
cone wavefunction which are "freed" during the collision, and this is indeed 
the picture we have in mind. However, we do not preclude that a significant 
number of these gluons may be created, in a light-cone gauge description, 
after the valence quarks of the colliding nuclei passed each other. In this 
case it is more proper to describe these gluons as produced by the collision. 

It would be nice to have a reliable calculation of the production cross 
section for freeing gluons in a central ion-ion collision. Unfortunately, this 
has not been achieved, up to now, even in the McLerran-Venugopalan model. 
There is an interesting numerical calculation of some of the features of early 
time production^, but so far it is not clear how to convert these calculations 
into an evaluation of the number of gluons produced at time f In order to 
proceed further we simply suppose that all the gluons in the light-cone wave- 
function having £^ ^ Q^ are freed. Thus we take as our initial distribution 
of freed gluons 



dN , dxG 

— c 



^Ql = c^i-^ Ql (16) 



d%dy d%d^ 

where c is a constant expected to be of order one. c cannot be calculated 
without knowing exactly which gluons are freed during the collision. At early 
times we can use free streaming kinematics 

z = t tank y (17) 



to get 



dN _ (iV^^-l)Qg t 



d%dz An'^aNc f^ - z^ 

Later on we shall argue that (18) should be a reasonable approximation for 
all times well before equilibration takes place. 

3.2 The dominant two-body scatterings at early times 

Immediately after the gluons in a given rapidity region are freed they begin 
to move freely along straight line trajectories until scattering with other 



freed gluons makes them deviate from their original direction. We suppose 
that the main scattering which occurs is gluon-gluon elastic scattering. Such 
an assumption should be justified for sufficiently large Qs- Consider those 
gluons having zero rapidity in a given frame of reference, where we suppose 
the z-direction is the colinear direction of the colliding ions. If gluons £i 
and £2 scatter into gluons d,'^ and d^ the kinematics of the scattering can be 
represented as 

4^ = (^10, ^ix, ^ly, diz) = iiil, 1, 0, 0) (19) 

4/. = ^2(1, -cosx, -sinx, 0) (20) 

i'l^ = i'i(l, cos9, sinOcosip, sinOsimp) (21) 

4, = (4 + ^2 - e,), (22) 

with i.'i determined by (^2)^ = 0. The initial gluons lie in the x,y plane and 
approach each other at an angle x with the initial gluon li being scattered 
by an angle 9 from its original direction into the final gluon d.'^. ip is an 
azimuthal angle, but here about the x-axis. Useful invariants are 

s={ii + hf = 2iMl + cos x) (23) 

- i = -{h - e[)^ = 2£ii[{ 1 - cos 6) (24) 

while (4 + £2-^1)^ = Ogives 



ni 



£1^2 (1 + COS x) 



(25) 



£1(1 — cos 9) + £2(1 + cosxcosQ + sin x sin 9 cos ip) 
We shall be dealing with small angle scattering {9 small) so that 

4 = ^'^'^^ + '"' ^^ 2 • (26) 

£2(1 + cos x) + 4^ sin X cos (p> + ^ (£1 - £2Cos x) 

In what follows we shall be following the change in transverse momentum 
of a gluon due to scattering with other gluons in the medium. For the gluon 

9 



initially labeled by £i the transverse momentum is £i_l = J ^Ix + ^\y = ^i- 
After scattering 

^i± = VW^)' + W.)' = ^i(l - y^^n'0) (27) 

in the small 9 approximation. Thus 

h± - ('i± = h- ([ + h— sin^cp (28) 

in the small 9 approximation. Using (26) it is straightforward to arrive at 

- 2£2sm\cosV + 4(1 + cosxfsin^4>\. (29) 

After averaging over the azimuthal angle according to 



^''-^/o 2^^^-^/o 2^<^--^-> ('"' 

one finds 

3.3 The rate of change of gluon transverse momentum 
at early times 

Now that we know how much transverse momentum a gluon loses during 
an early-time collision, given by (29) and (31), it is not difficult to give 
an expression for the rate of change of the transverse momentum of a zero 
rapidity gluon due to its initial small angle scatterings. The result is 

M^^ _ f dN .^ .da .^ ^^,^^ 



dt J d^xdH2' 'dn[ 

where 



/^foi-i^^-^^-''^''''^^^ P^' 



dn[ = sin9d9d^ (33) 
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I^i2| = V (^^1 - V2y - (t/lXt72)^- (34) 

Eq.32 expresses the fact that the rate of change of transverse momentum 
is given by the change of transverse momentum in a single scattering event 
times the rate of scatterings. The rate of scatterings is given by the cross 
section times the relative flux, {vul, times the density of scatterers. The 
density of scatterers is given by (18) as 

,N _ e(« - « ,N _ ,iAizile(« - il) (35) 



d^xdH2 ttQI d%dz An^aN^t 

where we approximate the £2— spectrum by a step function. | -1/121 can be 
evaluted from (19), (20) and (34) as 

I -1/12 1 = 1 + cosx (36) 

while, for small angle scattering, using (10) and (24), 
da aiVe 2 47r3 rf^2 ^ 

The dip in (37) and the lack of any (/? dependence in any of the other factors, 
beside A£i_|_, in (32) allow one to replace A£i_|_2f by A£i_|_ as given by (31). 
Thus 



(38) 



or 





dil^ aN^ rde^ rQs 

dt ~ ^ TTt J 9^ J 


di2{2h - 


dt 


3caNc 2 [^ de"^ 
27ct ^'Jei 6^ 


_caN 
nt 



1 d9^ 



(39) 



where the ^—integration is cut off at 9m, which quantity will be the topic of 
our next section. Finally, supressing the subscript 1 on ii±, one finds 

*' ^^"^^ Q? C„ 1/C. (40) 



dt 2'Kt 

a formula accurate in the logarithmic approximation 
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3.4 The minimuni scattering angle 

The cutoff on the 9— integration in (39) corresponds to a cutoff on the 
momentum transfer squared, t, given by 

— tm = H±^m = ^m (^l) 

with im the minimum momentum transfer allowed. One can view (39) as 
arising from the scattering of a particular hard gluon with the (relatively) 
soft gluon field due to all the other hard gluons which have been produced at 
a similar rapidity. However, this soft field should also be limited to size 1/g 
because of saturation. We can impose the saturation condition by requiring 
that our hard gluon have only one scattering with momentum transfer greater 
than or equal to i^ in a distance equal to l/im- Thus, im, and hence 9^ is 
determined by 

„—f^ J^ J AT 1 

(42) 



(43) 











/_ 


~*m da 
oo di 


dN 
d^x 


1 

f 


1. 


ing 


(10), 


(18) 


and 


(41), 


3ne finds 


















IT 


Ql 
















aN, 







TT 



(44) 



This is not quite the same screening mass suggested in Ref.l3. If one uses 
Eq.3 of Ref.l3 with ^^y given from (16) above and with all freed gluons 
taken in the integral of Eq.3 one finds 

When {aNcYQst is much bigger than one im, as given by (44), is larger than 
fiD, as given in (45). For our purposes the exact form of the cutoff in the 
infrared is not so important so long as the cutoff decreases no faster than 
some fractional power of Qgt. 
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3.5 The time dependence of £_i_ 

Now go back to (40). It is convenient to write this equation as 



_1_ dilit) 



dt 



caN.Q, ^^ iljt) 



Using (44), 

d{il{t)/Q 



aK, 



dinQst 



Tl 



nQst 



3 f TT 

in^ + inQst + M— rr 
2 Qj caNc 



Defining 



one finds 



^ = inQst, ^0= M-"""^" 



TT 



d{i±/Qs 
d^ 



aNr 



— c- 



TT 



[e - Co + Sin{iJQs)] 



(46) 



(47) 



(48) 



(49) 



which is our final equation. 

So far we have derived (49) when t, and hence ^, is not too large, and 
thus in the regime where i±/Qs is very close to 1. In the next section we 
shall argue that (49) is valid in a much wider range of times. Thus, we here 
give a brief discussion of the properties of the solution to (49) even when ^ 
corresponds to very large times. We are interested in finding the solution to 
(49) in the region ^ > with q±/Qs = 1 at .^ = 0. The dependence of i±/Qs 
on ^ is qualitatively clear. For C, not too large one may neglect the ini±/Qs 
term on the right-hand side of (49) to get 



{i±/Qs 



caNc 

271 



(e-Co 



eo 



(50) 



i±/Qs continues to decrease with increasing C,- When C, is large, and hence 
^±/Qs is small, the ^—dependence of i±/Qs is obtained by setting the right- 
hand side of (49) equal to zero. This gives 



{^±/Q 



sj ^ e 



-|«-?o) 



(51) 

The transition from (50) to (51) takes place when ^ is extremely close to ^i 
with ^1 given by 
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V caNc 
or, equivalently, when t is near ti given by 



ii = \H^ + a + Co (52) 



U = A^^^'\ (53) 

ttQs 

Equation (49) represents scatterings which transform transverse momenta 
into longitudinal momenta. At early times there is a negligible transfer of 
longitudinal momenta into transverse momenta. As we shall see in the next 
section it is exactly at the time when longitudinally moving particles have 
a significant probability of scattering and creating additional transverse mo- 
menta that equilibration is beginning to set in. 

There is another way to estimate when equilibration begins to set in. An 
equilibrated expanding gas in QCD should obey the ideal gas expansion law 

with T the (large) temperature of the gas. For an equilibrated system trans- 
verse momentum and temperature are proportional so one should have 

^^M . l^^jQ^y, (55) 

At early times, ^ — Co small, it is clear from (49) that the transverse mo- 
mentum in our expanding gluon system is not decreasing anywhere near fast 
enough to be close to equilibration. However, because of the lowering of im 
as t increases the rate of decease of i±/Qs grows with time and we may es- 
timate the time when equilibration is likely to begin to happen by equating 
(49) and (55). (In the next section we shall see that this estimate is the 
same as that coming from requiring that the mean free path of longitudi- 
nally moving gluons be less than the longitudinal length of the medium.) 
One finds 

- liUQs)' = - c^ [e - Co + Sin iJQs] . (56) 
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Taking {i±/Qs)'^ on the left-hand side of (56) from (50) and (temporarily) 
neglecting in i±/Qs one gets 

giving .^ = .^2 with ^2 determined by 



6 = J^ + Co - 3/2 = Ci - 3/2 (58) 

V cal\r 



and 



(^x/Q.)i = Sf-^. (59) 

One can now go back to (56) and check that neglect of the in i±/Qs term on 
the right-hand side of that equation amounts only to a very small correction 

to6- 

Thus, our guess is that equilibration starts to set in at a time ^2 given by 
(48) and (58) and at that time {i±/Qs)^ is given by (59). Also at ^2, from 

(44), 

iniQl/A')-iniil/A') = li^-^o). (60) 



Using (58) one sees that at ,^ = ,^2 




1 1 




aiQl) a{il) 



(61) 

so that running coupling effects should not be important in the analysis we 
have given. 

Finally, we note that 



1 . caA^c 
t2 = -pr-e'^ = —-—exp 



2-n 3 



caN, 2 



(62) 



Using the fact that parametrically Q? oc ^i\ R with R the nuclear radius, 
equilibration can take place before the one-dimensional expansion changes to 
a three-dimensional expansion. 
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3.6 From early times to equalibration 

Suppose we take as our initial reference frame the center of mass of the 
central ion-ion collision with the right-moving valence quarks crossing the 
left-moving valence quarks at t = 0. At ^ = gluons at zero rapidity are freed 
and begin to interact. At early times a gluon at zero rapidity can only interact 
with other gluons having small rapidity because all other gluons separate 
quickly in the z-direction from zero rapidity gluons. A gluon having rapidity 
y has a longitudinal velocity Vz = tanhy and can interact only with other 
gluons having rapidity close to y because it, too, separates longitudinally 
from all other gluons. Of course if two gluons having rapidity y interact we 
can change to a coordinate system, boosted by rapidity —y from our original 
system, where these gluons have zero rapidity and our calculation given above 
applies. Thus, our calculation for changes in transverse momentum, which 
is boost invariant, is also true for gluons having rapidity y. In particular, 
(49) remains valid for the rate of change of transverse momentum of a gluon 
having y if one simply identifies ^ as inQst—inchy and ,^o as £n{^^^) — inchy 
with (49) being valid for .^ > 0. 

To see when the approximation of scattering only between gluons having 
identical rapidities breaks down it is convenient to introduce ( by 



Then (49) and (50) become 

d{i±/Qsf caN, 



dC 



TT 



and 



ii±/Qs 



I 27r 
caNc 



I 27r 
caNr 



IcaNc 
2n 



+a-c. 



C + Mn{£±/Qs 



2tt ^ 



(63) 



(64) 



(65) 



with 



caNc 

TlQs 



27r 



tie 



(66) 
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So long as ''^^ ,q °. is small the approximation of having only gluons of equal 
rapidity interact is good. This is perhaps most easily seen by imagining that 
at time t a particular gluon begins to cross, making a finite angle with the z- 
axis, the gluon medium of thickness /\z ^ 2t. If our particular gluon is right- 
moving, then during a time t it will cross all the left-moving gluons in the 
medium. If during that time, corresponding to A(^ ^ 1, the gluon does not 
suffer enough interactions to change its momentum by a significant amount, 
then its important interactions must occur much later and with gluons having 
a similar rapidity since only such gluons have the same longitudinal velocity 
and so do not separate from our particular gluon which is essentially moving 
along a straight line trajectory. Now, using (64) and (65), one finds 



iSJQs) _ \ iiiJQs? _ 1 



(67) 



ii±/Qs) ii±/QsV 2C 

Thus so long as ( is large our procedure of calculation should be reasonable. 
However, somewhere when ( becomes on the order of 1 our whole procedure 
begins to break down as scattering between unequal rapdity gluons becomes 
important. In addition, at this time, scatterings which transform longitudinal 
momentum into transverse momentum are also becoming important, and 
they are not included in our formalism. 

Thus when ( is of order 1 we expect the interactions important for equi- 
libration to begin and our approximation to break down. As can be seen by 
comparing (58) and (63) this is essentially the same criterion we found in 
the last section. Thus, in our simple approach we are not able to follow glu- 
ons when the interactions essential for equilibration begin. We are, however, 
able to follow the system, on average, up to the time at which equilibrating 
interactions begin to occur. We have a good estimate of that time which 
is close to ti, given in (53), and of the transverse momentum of the gluons, 
given by (59), at the time when equilibration starts. However, because the 
density of gluons, given by (18), is of size c ^^l ~jv Ql^^^^ when equilibra- 
tion starts we expect that when equilibration is complete and ^ ~ ^j_ that 
{i±/Qsy ~ e~3 5i. It is encouraging that equilibration starts well before the 
expansion changes from one-dimensional to three-dimensional. 
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